#include"functions.h"

using namespace std;

int main()
{
  for(int n = 50; n < 101; n++)
    {
      double** A = new double*[n];
      for(int i = 0; i < n; i++)
	{
	  A[i] = new double[n]();
	  A[i][i] = 4.0;
	}
      for(int i = 0; i< n-1; i++)
	A[i][i+1] = A[i+1][i] = 1.0;
      Jacobi(A, pow(10, -10), n);
      std::cout << std::endl;
    }
    
  return 0;
}
